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Using the Two-dimensional Oscillation Program (TOP), we have explored the effects of rapid rotation on gravity modes 
in polytropic stars. Coriolis force, centrifugal distortion as well as compressible effects have been taken into account. 
Thanks to our complete calculation, we have first studied the validity domain of perturbative methods and started to 
explore properties of these modes. We focus on I = 1 in this analysis. 
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1. Introduction 

Rapidly rotating stars are common objects on the main 
sequence. These stars can be strongly distorted by the cen- 
trifugal force, and some of them are known to spin almost 
as fast as their break-up frequency. Fast rotators are easily 
found in some classes of pulsating stars, such as 8 Scuti, 
7 Doradus, f3 Cepheid or SPB (slow pulsating B) stars. 
Acoustic (p) modes are excited in 5 Set and [3 Cep stars, 
whereas gravity (g) modes are responsible for SPB and 7 
Dor oscillations. Rotation effects cannot generally be ne- 
glected for frequency computation even for slow rotation. 
However, when the rotation frequency f2 is small relative to 

mode frequencies w and to Q,k = 



GM/Rl q , rotation ef- 
fects can be treated as a perturbation (M and R eq are the 
mass and equatorial radius of the star). Such approaches 
ease frequency computations, but have a limited applica - 
tion range in term of rotation rate. Lignieres et al. (20061, 
Reese et al. (20061 have investigated the limits of pertur- 
bation theory for p-mode computation showing that taking 
the centrifugal distortion into account is needed to correctly 
predict mode frequencies. 

In this paper we investigate the effects of rapid rota- 
tion on gravity mode frequencies by performing a full two- 
dimensional computation of modes in polytropic models of 
stars. The effects of rotation on the modes are fully taken 
into account. 



2. Methods 

We consider for this work polytropic models of stars with a 
polytropic index N — 3 corresponding to a radiative zone. 
The 2D struct ure of the rotating po lytrope is computed 
as described in |Lignieres et al. ( 2006 1 . We have considered 
models spinning uniformly at the rotation frequency f2 be- 
tween and Q.7Q.K- Models are descri bed in a spheroida l 



This resolution is sufficient enough to accurately model the 
centrifugal effects for the maximal value of f2 we have con- 
sidered. Figure [l] shows the profile of the Brunt- Vaisala fre- 
quency N Q normalized by fl p K = GM / R^ ol (where R po i 

is the polar radius) for the two extreme models (O=0 and 
0.7f2if ).We see that N Q does not change too much in terms 
of values and profiles for the different models. 




Fig. 1. Solid black line shows the profile of Brunt-Vaisala 
frequencies N Q normalized by £l p K for the non-rotating star. 
Red lines show N a for the model il = Q.7Q.K along the polar 
(solid line) and equatorial radii (dashes). 



We solve the eigenvalue problem corresponding to adia- 
batic small perturbations of the polytropic model of star. To 
do that, we used t he Two-dimensiona l Oscillation Program 
(TOP) based 



on 



Reese et al. (2006) 



geometry (more details for instance in Reese et al. 2006). decomposed on spherical harmonics 



We have used a resolution of n r = 96 in the pseudo-radial 
direction and L max = 32 spherical harmonics in latitude. 



The modes being 
, the resolution is 
given by the truncation of the decomposition. For this work, 
we have considered spherical harmonics up to L = 40+ |m|, 
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i.e. 20 coupled harmonics. The radial discretization is the 
same as for the models, i.e. n r = 96. 

3. Results 

We have computed the solutions with I = 1 (m=— 1,0,1) 
and radial orders n = — 1... — 14 without rotation, and fol- 
lowed the evolution of the frequencies of these 42 modes 
by slowly increasing the rotation rate to 0.7£Ik- These fre- 
quencies are plotted in Fig. [2] 

From this data, we can compute the frequencies that 
would be determined in a perturbative approach. In such 
an approximation, frequencies are expressed as 

- pert 
UJ i — 
n.l, 



w„,i + c, m Q + cfn 2 



cl m Q 3 + o(n% (i) 



where uj — uj/SI p k and Cl — Q,/VL V K . We have considered 
approximations to the 3rd order. The term ct>„,z is noth- 
ing else but the frequency without rotation, whereas the 
coefficients C\ m are numerical computations of the ith- 
derivative d l uj n ,^ rn / 'd£l l at Q = obtained thanks to our 
full computation. Since these perturbative coefficients con- 
tain all the considered physics, they are ideal coefficients. 
In a sense, any perturbative theory aims to approximate 
the coefficients we have found. We have verified that the 
frequencies without rotation ui n / a re in perfect agreement 
(10 



& Mullan|l994 


1 . We have also 


ical expression 
within relative 


(Ledoux|1951l 

errors of 10 -9 



Concerning the coefficients 
Cf m and Cf m we have estimated the errors to be around 
10~ 5 . This last error is deduced from the errors estimated 
for our numerical derivatives. 

We have normalized frequencies by il p K since the polar 
radius is expected to be a slowly varying function of £1 in 
real stars, as opposed to R eq . 



error we allow on the frequency and ii) on the frequency 
£l p K that depends on the star. Figured shows the validity 
domain of 1st, 2nd and 3nd order methods for a typical 
7 Dor star (AT = 1.55Af0, R po i — I.QRq) for an error 
of 0.1 /iHz. In this case, perturbative approaches can be 
sufficient to determined the higher frequencies as long as 
£1 < 0.350^ but only for modes with uj > 2fL The do- 
main of validity would be more extended for larger error 
tolerances or for more massive stars (i.e. higher Sl p K ) such 
as SPB stars. Nevertheless we observe in any case a strong 
departure from the exact frequencies as uj < 20. This can 
be interpreted as the fact that the volume of the resonant 
cavity is modified by th e effects of the Coriolis force when 
uj < 2fl (see for instance Dintrans & Rieutord 2000). 
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Fig. 3. Period spacing AP,^ for I = 1 and m = as a 
function of the mode periods for different rotation speeds. 
Red crosses indicates sub-inertial modes. Dots indicate the 



asymptotic value from |Tassoul (19801. Periods are in unit 
of 2tt/Q^. 




Fig. 2. Evolution of the frequencies of I = 1 modes. 
Frequencies are computed in the co-rotating frame. 
Perturbative approximations have been tested for a star 
with M — 1.5Mq and R po i — 1.6Rq and an error 
of 0.1 ^Hz: green/red/blue parts of curves indicate that 
lst/2nd/3rd order is sufficient to reproduce complete cal- 
culations. The magenta line indicates uj = 2f2. 



For a non-rotating star, the period spacing AP nj ; = 
P n +i,l,o - P n ,l,o, where the period 

known to be -asymptotically- inde pendent of n, in the ab- 
sence of strong structure gradients ( |Tassoul|1980 1. Figure [3] 
shows the evolution of this spacing, for I = l,m = modes 
when rotation increases. While this spacing is almost con- 
stant at low rotation, even for low-order g modes, as soon 
as the rotation becomes non negligible, it is no more the 
case and the period spacing follows a more complex law. 
As a consequence, at higher rotation rates, Tassoul's for- 
mula cannot be used for mode identification. 

4. Conclusions 

This short report shows some preliminary results obtained 
on g modes. This work has been extended to I = 2 and 3 
modes and will be presented in more details in an forth- 
coming paper (Ballot ct al., in preparation). 
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We have then derived the validity domains for 1st, 2nd 
and 3nd order methods. These domains depend i) on the 
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